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Abstract 

We show that and how the Coulomb potential V(x) = Ze 2 /x can be regularized 
and solved exactly at the imaginary coupling Z e 2 . The new spectrum of energies is 
real and bounded as expected, but its explicit form proves totally different from the 
usual real- coupling case. 



1 Introduction 



Quantum mechanics often works with the exactly solvable simplified models. For the 
precise fits of data or for some more subtle quantitative analyses, unfortunately, the 
number of solvable models is too limited. In D dimensions, the only really useful and 
easily tractable interactions are harmonic oscillators and/or the central Coulomb well 
V^ z \r) = —Ze 2 /\r\. A new way out of this deadlock emerges within the framework 
of the alternative, so called "VT symmetric" quantum mechanics. With its complex 
Hamiltonians H breaking both the parity V and the time-reflection symmetry T 
and commuting only with their product VT, this formalism was proposed by Bessis 
|I[ and by Bender et al 0, 0] as a possible way towards weakening of the standard 
requirements of Hermiticity. 

Several new exactly solvable VT symmetric models have been proposed recently 
Q. This is a promising development with possible applications ranging from field 
theories || to supersymmetric models || and from quasi-classical methods [^J to 
perturbation theory [[/]. 

Even the solvable harmonic oscillator itself acquires a richer spectrum after its 
consequent VT symmetric regularization in D dimensions [§]. The detailed struc- 
ture of spectrum of this prominent example does not in fact offer any really serious 
surprise. A manifest violation of the parity V is compensated by an emergence of the 
so called quasi-parity q — ±1 tractable as a signature of two equidistant subspectra. 
The new quantum number q degenerates back to the eigenvalue of parity after a 
return to the standard Hermitean and one-dimensional oscillator. 

No immediate surprise emerges also for the quartic anharmonic oscillator ||. 
The situation only becomes less clear after one moves towards the asymptotically 
vanishing models. They exhibit several counterintuitive properties and open new 
mathematical challenges pTD | . In particular, the popular Coulomb potential did not 
not even seem particularly suitable for any immediate VT symmetric regularization 



1 



11]. A psychological barrier has been created by the numerical and semiclassical 
studies of the general power-law forces V(x) ~ — (ix) s . They may be well defined 
everywhere near the harmonic exponents 5 = 2, 5 = 6 etc fl2fl . At the same time, the 
related analyses hinted that it is apparently difficult to move beyond the Herbst's 
singularity located at 5 = 1 JT3|. 

In what follows, we intend to employ a slightly different strategy and try to study 
the Coulomb problem directly, via its well known correspondence to the harmonic 
oscillator. This correspondence is based on an elementary change of variables. Its 
background dates back to the nineteenth century mathematics and, in particular, to 



the work of Liouville | Hf . The Newton's monograph [|15| cites also Fivel Jl6| as a 
newer source of the idea. In the contemporary literature (cf., e.g., [^] for further 
references) people usually speak about the Kustaanheimo - Steifel (KS) transfor- 
mation fl8| . In all the implementations of this idea the parameters appearing in 
the Coulombic and oscillator problems are interrelated of course. Details will be 
mentioned below. Preliminarily, let us only warn the reader that all the KS-type 
mappings can also change the dimensions and angular momenta and that the en- 
ergies of one problem are related to the coupling constants of the other one and 
vice versa. Within the "normal" quantum mechanics, all this has already been thor- 
oughly discussed elsewhere: In ref. for D = 3 and in ref. |2(J for the continuous 
transformation between Coulomb problems and harmonic oscillators in various di- 
mensions. 



2 Liouvillean changes of variables 



The change-of-variable approach to the Coulombic bound-state problem enables us 
to start directly from the harmonic oscillator potential W(r) = r 2 or, in the present 
less traditional context, from its VT symmetric radial Schrodinger equation 

dr 2 



x (r) = e 2 x{r) 



(1) 
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of ref. |J, using the complex coordinate r = x — ic with real x G (— oo, oo) and with, 
say, positive c > 0. This means that the integration path has been shifted down 
from the position where it would cross the strong centrifugal singularity. Such a 
regularization preserves the asymptotic decrease of the normalizable solutions. Only 
in the limit c — > and beyond the trivial one-dimensional case one has to omit 
all the so called irregular solutions |ZT| defined by their "physically unnacceptable" 
x{r) ~ r~ l behaviour near the origin. 

Within the framework of the general Liouville method the change of variables 
mediates a transition to the different potential V(t). It is easy to show that once 
we forget about boundary conditions one merely has to demand the existence of an 
invertible function r = r(t) and its few derivatives r'(t), r"(t), .... Then, the explicit 
correspondence between the two bound state problems may be explicitly given by the 
elementary formulae. From our original eq. (0) (i.e., in our case, harmonic oscillator) 
one obtains the new (i.e., in our case, Coulombic) Schrodinger equation 



dt 2 



with the new wave functions 



t 2 



(2) 



*(*) = xW)]/\p(t) 



(3) 



and with the new interaction and the new energies [[14 
L(L + 1 



t 2 



+ V(t)-E=[r\t)] 2 i^^ + W[r(t)]-^ 



H 


>"(*)" 


2 1 


>'"(*)" 


r'{t) 


2 


_ r'(t) _ 



Thus, it only remains for us to re-analyse the boundary conditions. 



3 VT symmetric KS transformation 

Without any serious formal difficulties let us extend the scope of the present consider- 
ations to all the singular forces W(r) = W{r) + / jr 2 and and/or V(t) = V(t) + F/t 2 . 



3 



Both these central forces may act in the respective d and D dimensions. This means 
that 

1(1 + l) = {[ J + (d- 3)/2] [j + (d- l)/2] + /} 

or rather 

(/ + 1/2) 2 (= a 2 ) = [j + (d-2)/2}] 2 + f 

(with Re a > 0) and 

(L + 1/2) 2 (= A 2 ) = [J +(D — 2)/2]] 2 + F 

(with Re A > 0) where the partial waves are numbered by the respective integers 
j = 0, 1, . . . and J = 0, 1, . . .. 

An important simplification of our effort is provided by our knowledge of the 
complete harmonic oscillator solution as derived in ref. ||. Its two equidistant 
subsets of energies 

£<2 — e \n, q ) = 4™ + 2 - 2qa, q = ±1, n = 0, 1, . . . 

correspond to the two families of the Laguerre-polynomial wave functions 

Xin, q) {r)=Nr l l 2 -i«e- r2 ' 2 L^\r 2 ). 

Integration path r = x — ic lies in the lower half of the complex plane and does 
not change after the subsequent V and T transformations r = r[x) — * — r and 
— r — > —r* = r(—x). 

In the spirit of the above-mentioned KS mapping of harmonic oscillators on 
Coulombic bound states we now have to define a complex variable t as a re-scaled 
square of r(x) such that the resulting path t(x) remains VT invariant. Our require- 
ment implies that the complex plane of r will cover twice the complex plane of t. 
In such an arrangement, our lower half plane of r should cover the Riemann sheet 
given as a whole plane of t which is cut upwards from the origin. In the polar rep- 
resentation, one has r ~ exp(— i<p) mapped upon t ~ exp(— 2icp) with 6 (0, n). 
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Once we introduce a suitable free parameter k > we can put, say, r 2 = 2k 2 z and 
then rotate the z— plane (which is cut, by construction, along the real and positive 
semi-axis) by the angle 7r/2 giving t = i z. Our final recipe 



r z = -2in z t (4) 

maps the above-mentioned straight line r(x) = x — ic upon a curve t(x) = u + iv. Its 
real part u = u(x) = xc/n 2 and imaginary part v = v(x) = (x 2 — c 2 )/2n 2 form, as 
required, a VT symmetric and upwards-oriented parabola v = — c 2 /2 k 2 + (k 2 /2c 2 ) u 2 
in complex plane. Obviously, a small asymptotic deformation of our original curve 
r(x) with the modified shifts c = c(x) ~ l/x l+v would transform the parabola to a 
pair of lines which are parallel to the cut in the asymptotic domain of \x\ ^> 1. 

Having achieved a VT symmetry in the complex plane of t, we may move to the 
(trivial) insertions and conclude that all the above-mentioned harmonic oscillator 
bound-state solutions are in a one-to-one correspondence to the solutions of the 
Coulombic Schrodinger equation (0), 



d 2 L(L + 1) Ze 2 



= EW(t) , t = u{x)+iv(x), xeM.(B) 



dt 2 t 2 t 
The underlying assignment of constants is such that a = 2 A while k itself becomes 
n— dependent, k 2 = 2 Z e 2 /e 2 = Z e 2 / (2n + l — 2q A). In full detail one gets the new, 
Laguerre-polynomial wave functions 

*(«,,)(*) = Mt l ' 2 - qA e iK2t L^ 2(lA \-2iK 2 t) (6) 

and their energy spectrum specified by the elementary formula 

E ^ = K U)= (2n + f_ e 2gA)2 9 = ±1 . « = 0,1,.... (7) 
This is our main result. 



4 Discussion 
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4.1 Consequences of the curvature of our integration path 

The latter two formulae exhibit several unusual features. The first concerns the 
asymptotics of the wave functions which are determined by the decreasing exponen- 
tial exp(in 2 t). Its form re-confirms the correctness of the above, slightly counter- 
intuitive KS-dictated choice of our VT symmetric integration path. Asymptotically, 
it encircles more or less closely the positive imaginary axis in t plane. This clarifies 
the apparent paradox. 

The second unexpected result is the positivity and unusual n— dependence of the 
energies. This can be related to the choice of the KS integration path t(x) again. In 
the very vicinity of the origin, one can visualize this path as a circle with radius a, 

u 2 (x) + v 2 (x) = a 2 , \x\ < 1. 

From the appropriate definitions we get the formula 

a = c 2 (0)/2n 2 M + O(x 2 ) 

and see that this radius is n— dependent and increases with the growth of this prin- 
cipal quantum number, a ~ nc 2 (0)/Ze 2 , n ^> 1. As a consequence, an "effective 
charge" of our VT symmetric Coulomb potential appears to decrease with n since 

— + 0(t) = 0(l/n). 
a 

This offers a "rule-of-thumb" guide to the unusual and certainly counterintuitive 
n— dependence of the energy levels (|7[). Of course, in practice, a preferred integration 
path will be n— independent. In such a case, the n— dependence re-appears in the 
small- x deformation of the initial harmonic-oscillator path with c = c n (x) = 0(l/n). 
Such a flexible transfer of the excitation-dependence throws also a new light on the 
complexified KS transformation itself. 



iZe 2 
t 
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4.2 "Flown-away" energies and unavoided crossings 

Let us return in more detail to the A— dependence of our energies (0). Firstly, we 
notice their power-law dependence on n and A (with exponent = -2) as somewhat 
similar to the spectra in VT symmetric oscillator well || (with exponent = +1) and 
in the Morse potential of ref. |22| (with exponent = +2). 

In the present case, obviously, we have to distinguish between the two separate 
families -E( n , g ) with q = +1 (cf. Figure 1) and q — — 1 (cf. Figure 2). The latter 
set is, up to its sign, analogous to the ordinary Coulombic spectrum. By far not so 
the former one. Its energies enrich and dominate the spectrum. The n^—th energy 
"flies away" and disappears from the spectrum at A div = n div + 1/2. Moreover, at all 
the positive integers and half-integers A one encounters the unavoided level crossings. 
In contrast to the harmonic case, they appear at both the opposite and identical (viz., 
positive) quasiparities q. The former case takes place at A = A cr u = (n — n')/2 while 
in the latter case we must fulfill the condition A = A cr a = (n + n' + l)/2. A sample 
of this phenomenon is given in Figure 3. 

Formally the unavoided crossings generate certain identities which connect dif- 
ferent Laguerre polynomials (cf. their sample in ref. ||). In applications, these 
"critical" cases are not exceptional at all. For F = forces without a spike, the 
critical integer or half-integer coordinates A crit = J — 1 + D/2 correspond precisely 
to the physical (namely, integer) dimensions D and partial waves J. 

In the conclusion, let us not forget about many open questions. Pars pro toto, 
one could mention a not yet clear possibility of re-interpretation of our bound states, 
say, in the limit t — > real, i.e., beyond the mathematical and apparently natural 
boundaries of our present approach. Moreover, we must keep in mind that via our 
complexification of the coordinates we of course broke their immediate connection 
to any standard D— dimensional problem. In this sense, our F = and F / 
Hamiltonians differ just in an inessential way. One could even prefer the latter, 
Kratzer-like option as a model which is formally simpler, due to the generic absence 
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of the puzzling unavoided crossings. At F ^ the structure of the spectrum of our 
present Coulomb model becomes also richer and, in this sense, more interesting. 
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Figure 1. Quasi-even levels (n,q) (units Z=e=l) 




Figure 2. Quasi-odd levels (n=0-5, units Z=e=l) 
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Figure 3. Sample of unavoided crossings 




